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Further X-ray observations of EXO 0748—676 in 
quiescence: evidence for a cooling neutron star crust 
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ABSTRACT 

In late 2008, the quasi-persistent neutron star X-ray transient and eclipsing binary 
EXO 0748-676 started a transition from outburst to quiescence, after it had been 
actively accreting for more than 24 years. In a previous work, we discussed Chan- 
dra and Swift observations obtained during the first five months after this tran- 
sition. Here, we report on further X-ray observations of EXO 0748-676, extending 
the quiescent monitoring to 1.6 years. Chandra and XMM-Newton data reveal qui- 
escent X-ray spectra composed of a soft, thermal component that is well- fitted by a 
neutron star atmosphere model. An additional hard powerlaw tail is detected that 
changes non-monotonically over time, contributing between 4 and 20 percent to the 
total unabsorbed 0.5-10 keV flux. The combined set of Chandra, XMM-Newton and 
Swift data reveals that the thermal bolometric luminosity fades from ~lx 10 34 to 
6 x 10 33 (D/7.4 kpc) 2 erg s _1 , whereas the inferred neutron star effective tempera- 
ture decreases from ~ 124 to 109 eV. We interpret the observed decay as cooling of 
the neutron star crust and show that the fractional quiescent temperature change of 
EXO 0748-676 is markedly smaller than observed for three other neutron star X-ray 
binaries that underwent prolonged accretion outbursts. 
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1 INTRODUCTION 

EXO 0748-676 is an intensively studied low-mass X-ray bi- 
nary that was initially discovered with the European X- 
ray Observatory SA Tellite (EXOSAT) in 1985 February 
IjParmar et alj|l985f ). However, in retrospect the source al- 
ready appeared active in EXOSAT sl ew survey observation s 
several times beginning 1984 July IjRevnolds et al.1 1 19991) ■ 
whereas the earliest detection dates back to 1980 May, 
when EXO 0748-676 was serendipitously observed with the 
EINSTEIN satellite (|Parmar et alj|l986h . The system ex- 
hibits irregular X-ray dips and displays eclipses that last for 
~ 8.3 min and recur every 3.82 hr, which allow the unam- 
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biguous determinat i on of the orbital period of the binary 
(|Parmar et al.lll986l : IWolff et al.ll2009h . 



The dete c tion of type-I X-ray bursts (e.g., 
iGottwald et al.1 Il986f ) conclusively identify the com- 
pact primary as a neutron star. A few X-ray bursts have 
been observed that exhibited photospheric radius expansion 
(PRE), which indicates that the Eddington luminosity is 
reached near the burst pe ak and allows for a distance esti- 
mate towards the source (| Wolff et al.ll2005l ; iGallowav et al.1 
120081 ). For a Helium-dominated photosphere, a distance of 
D — 7.4 ±0.9 kpc can be derived, while assuming solar com- 
position results in a d istance estimate of D = 5.9 ± 0.9 kpc 
IjGallowav et al.ll2008r ). The rise time and duration of the 
PRE bursts observed from EXO 0748-676 suggest pure 
Heliu m ignition, rendering 7 .4 kpc as the best distance esti- 
mate (|Gallowav et al.ll2008l ). although this value is subject 



2 N. Degenaar et al. 



37 (D/7.4 kpc) 2 erg s 1 by various satellites for over 
Similar prolonged accretion episodes continuing 



to sev eral uncertainties l|Wolff et al.l 120051 ; iGallowav et al.l 
l2008bl ). 

At the time of its discovery, EXO 0748-676 was 
detected at 2-10 keV lu minosities of ~ (1 — 7) x 
10 36 (D/7.4 kpc) 2 ergs" 1 (|Parmar et al.lll986l ). However, 
during the EINSTEIN observation of 1980, several years 
prior to the EXOSAT detections, it displayed a 0.5-10 keV 
luminosity of ~ 5 x 10 33 (D/7.4 kpc) 2 erg s" 1 (jParmar et all 
1 19861 ; iGarcia fe Callananl Il999l ) . The source can therefore 
be classified as a transient X-ray binary. Nevertheless, 
such systems typically exhibit accretion outburst s that last 
only weeks to months (e.g., I Chen et alj Il997l ). whereas 
EXO 0748-676 was persistently detected at luminosities of 
~ 10 
24 years 

for years to decades have been observed for a few other sys- 
tems, which are termed quasi-persistent X-ray binaries (e.g., 
IWiinandslliool ). 

In 2008 August-September, observations with the Pro- 
portional Counter Array (PCA) onboard the Rossi X- 
ray Timing Explorer (RXTE) and Swift's X-ray Telescope 
(XRT) in dicated that the X- ray flux of EXO 0748-676 was 
declining (|Wolff et al.ll2008al lbh. Optical and near-IR obser- 
vations of the optical counterpart, UY Vol, performed in 
2008 October showed that the optical em ission had also 
faded compared to the b righter X-ray state (|Hvnes fc Jones! 
|200S| , |2009| ; iTorres et J] 120081 ) ■ These events indicated that 
the accretion was ceasing and that the system was transi- 
tioning from outburst to quiescence. This is also illustrated 
by Fig. [1] which displays the X-ray lightcurve of EXO 0748- 
676 as observed with the All-Sky Monitor (ASM) onboard 
RXTE since 1996. The decrease in source activity is clearly 
seen around ~ 4600 days. 

Chandra observations carried out in 2008 mid-October 
(i.e., after the transition to quiescence started) revealed an 
X-ray spectrum composed of a soft, thermal component 
joined by a hard power law tail that domina tes the spec- 
trum above ^2 — 3 keV (JDegenaar et al.ll2009l . see also Sec- 
tion [2]4]). This is freque ntly seen for neutron star X-ray bina- 
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nesin quiescence (e.g., Kutlcdgc et al. 1999; in t /.an d et al.l 
l200ll ; iTomsick et al.l 120041 ). The non-thermal component is 



usually well-fitte d by a simple powerlaw with index 1-2 (e.g., 
lAsai et al.1 Il996l ). The fractional contribution of the hard 
powerlaw tail to the 0.5-10 keV X-ray flux widely varies 
amongst sources and possibly als o with changing luminosity 
jjonker et al.ll2004l ; |Jonkerll2008l ). The physical process that 
is responsible fo r the powerlaw spectral component rem ains 
elusive (see e.g.. lCampana et al.lll998l : ICampanall2003h . 

Although the soft spectr al component has be en as- 
cribed to low- level accretion l|Zampieri et al.l |l995), it is 
most often interpreted a s thermal surface radiation from 
the cooling neutron star l|Brown et al.l 1 19981 ). According to 
this model, the accretion of matter compresses the neu- 
tron star crust, which induces a series of electron captures, 
neutron emissions and pycnonuclear fusion reactions (e.g., 
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lHaensel fc Zdunidl990l . l2003l . l2008l : lGupta et al.ll2007l ). The 
heat energy released in these processes is spread over the 
neutron star via thermal conduction. 

The neutron star cools primarily via neutrino emissions 
from the stellar core, as well as photon radiation from the 
surface. The former depends on the equation of state of cold 
nuclear matter and the central density of the neutron star 
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Figure 1. RXTE/ ASM 20-day averaged lightcurve (1.5-12 keV) 
of EXO 0748-676, illustrating the cessation of the outburst in 
2008 August/September. For reference: the dashed vertical line 
corresponds to 2008 December 31. The arrows indicate the times 
of our four sequences of Chandra observations, which were per- 
formed when the source dropped below the detection limit of 
RXTE (both of the ASM and the PCA). 



(e.g.. lYakovlev fc Pethick1l2004l ; IPaee et al.ll200rJ ). The neu- 
tron star core reaches a thermal steady state in ~ 10 4 years, 
yielding an incandescent emission from the neutron star sur- 
face set by the time-averaged accretion rate of the system, 
as well a s the rate of neutrino emissions fro m the stellar 
core fe.g. jBrown et al.ll 19981 ; IColpi et al.ll200ll ). When com- 
bined with estimates of the outburst history, observations 
of quiescent neutron stars can constrain the rate of neutrino 
emissions, thereby providing i nsight into the inte rior prop- 
erties of the neutron star fe.g.. lHeinke et al.ll2009i ). 

Once the steady state is reached, the neutron star core 
temperature will not change appreciably during a single out- 
burst, but the temperature of the crust can be dramatically 
altered. In regular transients that have a typical outburst 
duration of weeks to months, the crustal heating processes 
will only cause a s light increase in the crust temperature 
IjBrown et al.lll9 98). However, in quasi-persistent X-ray bi- 
naries the prolonged accretion episodes can cause a signif- 
icant temperature gradient between the neutron star crust 
and core. Once the accretion ceases, the crust is expected 
to thermally relax on a time scale of years, until equilib - 
rium with the core is re-established (|Rutledge et al.l 12002). 
During the initial stages of the quiescent phase the thermal 
emission will therefore be dominated by the cooling crust, 
whereas eventually a quiescent base l evel is reached that is 
set by the thermal s tate of the core IjWiinands et a.1.1 l200ll : 
iRutledge et al.1 -2002). This provides the special opportunity 
to separately probe the properties of the neutron st ar crust 
l|Haensel fc ZdunikH2008l ; iBrown fc Cummindl2009l ) . 

In 2001, the neutron star X-ray binaries KS 1731-260 
and MXB 1659-29 both made the transition to quiescence, 
following accretion episodes of 12.5 and 2.5 years, respec- 
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lively dWiinands et al.ll200ll , |2002| , 120031 . 12004 ICackett et all 
12009 , 120081 ). More recently, in 2007, the ~ 1.6-year long out- 
burst of XTE J1701-462 came to a halt ll Altamirano et ail 
120071 ; iHoman et alll2007l ; iFridriksson et alj|2010l ). All three 
systems were subsequently monitored with Chandra and 
XMM-Newton, which revealed that thermal flux and neu- 
tron star temperature were gradually decreasing over the 
course of years (see also Section [4]). This can be interpreted 
as cooling of the neutron star crust that has been heated 
during the prolonged accretion outburst. Successful mod- 
elling of the observed quiescent X-ray lightcurves with neu- 
tron star thermal evolution models supports this hypothe- 
sis and provides important constraints on the crust proper- 
ties, such as the thermal conductivity (|Shternin et al.ll2007l ; 
iBrown fc Cummindl2009l ). 

Along these lines we have pursued an observational 
campaign of EXO 0748-676 to study the time evolution 
of the quiesc ent X-ray emiss i on fol lowing its long accretion 
outburst. In iDegenaar et al.l 1)20091 ) . we discussed Chandra 
and Swift observations obtained between 2008 September 
28 and 2009 January 30. We found a relatively hot and 
luminous quiescent system with a temperature of kT^§ ~ 
0.11-0.13 keV and a thermal 0.01-100 keV luminosity of 
~ (8-16) x 10 33 (D/7.4 kpc) 2 erg s _1 . No clear decrease in 
effective temperature and thermal bolometric flux was found 
over the five-month time span. In this paper we report on 
continued Swift and Chandra observations of EXO 0748- 
676 during its quiescent state. In addition, we include an 
archival XMM-Newton observation performed ~ 2 months 
after the cessation of the outb urst. Previous Chandr a and 
Swift observations discussed bv lDegenaar et al.l l|2009l ) were 
re-analysed in this work in order to obtain a homogeneous 
quiescent lightcurve. 



2 OBSERVATIONS AND DATA ANALYSIS 

Table [TJ gives an overview of all new observations of 
EXO 0748-676 discussed in this paper. A list of earlier Chan- 
dra and Swift observ ations obtained during the quiescent 
phase can be found in lDegenaar et all (1 20091 ). 



2.1 XMM-Newton 

EXO 0748-676 was observed with the European Photon 
Imaging Camera (EPIC) onboard XMM-New ton on 2008 
Nove mber 6 from 08:30-16:42 ut (see also iBassa et al.l 
I2009J). The EPIC instr ument consists of tw o MOS detec- 
tors ([Turner et alj|200ll ) and one PN camera IjStruder et al.l 
l200ll ). which are sensitive in the 0.1-15 keV energy range 
and have effective areas of 922 cm 2 and 1227 cm 2 (at 1 
keV), respectively. Both the PN and the two MOS instru- 
ments were operated in full window mode and using the 
medium optical blocking filter. Data reduction and analysis 
was carried out with the Science Analysis Software (SAS; v. 
9.0.0). We reprocessed the Original Data Files (ODF) using 
the tasks emproc and epproc. To identify possible peri- 
ods of high particle background, we extracted high-energy 
lightcurves (> 10 keV for the MOS and between 10-12 keV 
for the PN). No strong background flares occurred during 
the observation. The net exposure times are 29.0 and 22.9 ks 



Table 1. Observation log 



Satellite 


Obs ID 


Date 


Exp. time 
(ks) 


XMM 


0560180701* 


2008-11-06 


29.0 (MOS) 
22.9 (PN) 


Swift 


31272016 


2009-02-13 


3.5 


Swift 


31272017 


2009-02-20 


1.1 


Swift 


31272018* 


2009-02-23 


5.1 


Chandra 


9071* 


2009-02-23/24 


15.8 




10871* 


2009-02-25 


9.6 


Swift 


31272019* 


2009-03-01 


3.2 


Swift 


31272020 


2009-03-10 


5.1 


Swift 


31272021 


2009-03-16 


4.6 


Swift 


31272022 


2009-04-09 


3.5 


Swift 


31272023 


2009-04-16 


2.8 


Swift 


31272024 


2009-04-23 


1.8 


Swift 


31272025 


2009-05-07 


4.5 


Swift 


31272026 


2009-05-14 


3.6 


Swift 


31272027 


2009-05-28 


3.4 


Swift 


31272028 


2009-06-05 


4.1 


Chandra 


9072* 


2009-06-10 


27.2 


Swift 


31272029* 


2009-06-11 


4.3 


Swift 


31272030 


2009-06-18 


3.9 


Swift 


31272031* 


2009-06-26 


5.5 


Swift 


31272032* 


2009-07-03 


4.8 


Swift 


31272033 


2009-07-18 


5.5 


Swift 


31272034* 


2009-07-25 


5.8 


Swift 


31272035* 


2009-07-31 


10.3 


Swift 


31272036* 


2009-08-15 


9.4 


Swift 


31272037 


2009-08-25 


1.1 


Swift 


31272038 


2009-08-26 


7.4 


Swift 


31272039* 


2009-09-08 


4.7 


Swift 


31272040* 


2009-09-09 


4.3 


Swift 


31272041 


2009-10-01 


1.9 


Swift 


31272042 


2009-10-02 


1.8 


Swift 


31272043 


2009-10-07 


2.0 


Swift 


31272044 


2009-10-08 


2.1 


Swift 


31272045 


2009-10-09 


2.3 


Swift 


31272046* 


2009-11-05 


-1.2 


Swift 


31272047* 


2009-12-21 


9.4 


Swift 


31272048* 


2010-10-01 


9.6 


Swift 


31272049 


2010-02-12/13 


11.3 


Swift 


31272050* 


2010-03-12/13 


9.5 


Chandra 


11059* 


2010-04-20 


27.4 



Note. - The observations marked with an asterisk contain 
(part of) eclipses. The listed exposure times represent the 
duration of the observations uncorrected for eclipses. 



for the MOS and PN, respectively. EXO 0748-676 is de- 
tected at count rates of 0.16 ± 0.01 counts s" 1 (MOS) and 
0.55 ± 0.01 counts s" 1 (PN). 

Source spectra and lightcurves were obtained with the 
software task evselect, using a 35 arcsec circular re- 
gion and applying pattern selections 0-12 and 0-4 for the 
MOS and PN data, respectively. Corresponding background 
events were extracted from a circular region with a radius 
of 70 arcsec. For the MOS cameras, the background was po- 
sitioned on a source-free region on the same CCD as the 
source. For the PN instrument, the background events were 
extracted from an adjacent CCD, at the same distance from 
the readout node to ensure similar low-energy noise. The an- 
cillary response files (arf) and redistribution matrices (rmf) 
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were generated for each of two MOS and the PN cameras 
with the tasks arfgen and rmfgen. 

The EPIC li ghtcurves show two full eclipses (see also 
iBassa et al]|2009l ) , corresponding t o eclipse cycles 54384 and 
54385 in the numbering system of IParmar et al.l (J1986T) . To 
calculate the correct non-eclipse time-averaged fluxes, we re- 
duce the exposure times for each instrument by 500 s per 
eclipse, which i s the approximat e length of the eclipses of 
EXO 0748-676 (| Wolff et al .1120091) Fl Using the tool grppha, 
the spectra were grouped to contain a minimum of 20 pho- 
tons per bin. 



2.2 Chandra 

We obtained three new Chandra observations of EXO 0748- 
676 using the S3 chip of the Advanced CCD Imaging Spec- 
trometer (ACIS; ICarmire et al.l 120031 ). The ACIS detector 
is sensitive in the 0.1-10 keV passband and has an effective 
area of 340 cm 2 at 1 keV. The first observation consists of 
two separate exposures obtained on 2009 February 23-24 
22:07-03:15 UT (obs ID 9071) and 2009 February 25 12:32- 
15:59 UT (obs ID 10871), lasting for ~ 15.8 and ~ 9.6 ks, 
respectively. In both data sets, EXO 0748-676 is clearly de- 
tected at a count rate of 0.17 ± 0.01 counts s . This is a 
factor ~ 1.5 lower than observed in 2008 October, when 
the source was detected with Chandra / ACIS-S at a rate 
of 0.24 ± 0.01 counts s _1 . Two full eclipses are seen in the 
lightcurve of observation 9071, while one eclipse is present 
in that of 10871 (eclipse cycle numbers 55071, 55072 and 
55080, respectively). 

A second Chandra observation was carried out on 2009 
June 10, from 12:36-21:16 UT, with an exposure time of 
27.2 ks (obs ID 9072). In this observation EXO 0748-676 is 
detected at a count rate of 0.16 ± 0.01 counts s _1 and the 
lightcurve shows two full eclipses (cycles 55740 and 55741). 
Furthermore, a 27.4 ks exposure was taken on 2010 April 20 
from 02:37-11:28 UT (obs ID 11059), which captured three 
full eclipses (see Fig. [2] these correspond to eclipse cycle 
numbers 57708, 57709 and 57710), and detected the source 
at a count rate of 0.14 ± 0.01 counts s _1 . Similar to our 
treatment of the XMM-Newton data, we reduce the expo- 
sure times of all Chandra observations by 500 s per eclipse. 
There are no indications of background flares, so the full 
data set was used in further analysis. 

We reduced the data employing the CIAO tools (v. 4.2) 
and standard Chandra analysis threads. For all three ob- 
servation sequences, the ACIS-S3 CCD was operated in a 
1/8 sub-array, resulting in a frame-time of 0.4 s. For the ob- 
served count rates, the pile-up fraction is < 2 percent, so 
no further corrections were made. Spectra were extracted 
with the tool psextract, while the rmf and arf files were 
created using mkacisrmf and mkarf, respectively. We em- 
ploy a circular region of 3 arcsec to obtain source events and 
a 10-25 arcsec annulus for the background. We also repro- 
cess ed the Chandra obser vation obtained in 2008 October 
(see iDegenaar et al.l 120091 ') to benefit from the calibration 
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1 As shown by I Wolff et al.l (2009), the duration of the eclipses of 
EXO 0748-676 varied between ~ 484 and 512 s over the years 
1996—2008. These small uncertainties in the eclipse duration do 
not affect our results. 
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Figure 2. Chandra/ ACIS-S 0.5-8 keV lightcurve of EXO 0748- 
676 obtained on 2010 April 20 (obs ID 11059). Each point repre- 
sents 200 s of data. Three eclipses are visible. 

update that was released in 2009 December. Prior to spec- 
tral fitting, the spectra were grouped to contain at least 20 
photons per bin. 



2.3 Swift 

In addition to the Chandra and XMM-Newton observations, 
we have been monitorin g EXO 0748-67 6 on a regular ba- 
sis with the XRT (|Burrows et all I2005J ) aboard the Swift 
satellite. The instrument has an effective area of 110 cm 2 at 
1.5 keV and is operated in the energy range of 0.2-10 keV. 
Starting in 2008 late-September, approximately 2 — 3 point- 
ings were performed each month with a typical duration of 
~ 3 — 5 ks per observation, and a separation of ~ 1 — 2 weeks. 
From 2009 November onwards, the cadence was lowered to 
one observation per month with a longer exposure time when 
possible (see Tabled]). EXO 0748-676 is detected in the XRT 
observations at count rates of ~ (1 — 5) x 10 -2 counts s _1 . 

All Swift/XRT observations were obtained in the 
photon-counting (pc) mode and were processed using the 
xrtpipeline with standard quality cuts (event grade 0-12). 
Using xselect (v. 2.4), we extracted source spectra from 
a circular region with a radius of 35 arcsec (~ 15 pixels), 
which opti mises the signal to noise ratio at the observed 
count rates l|Evans et al.l 120071 ). Corresponding background 
events were averaged over three source-free regions of sim- 
ilar shape and size. Employing the tool XRTEXPOMAP, we 
created exposure maps to account for the effective area of 
the CDD, while arfs generated with xrtmkarf account for 
vignetting and point-spread-function corrections. The latest 
rmf (v. 11) was obtained from the CALDB database. 

Due to low statistics, it is not possible to identify 
eclipses in t he Swift lightcurv es. Therefore, we used the 
ephemeris of I Wolff et al] (|2009l ) to determine during which 
observations eclipses were occurring (see Table [TJ . To cal- 
culate the correct non-eclipse time-averaged fluxes, the ex- 
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Figure 3. Spectra of the Chandra observations of 2008 October 
(black) and 2010 April (red), along with the model fits (solid 
lines). The separate contributions of the NSATMOS and POWERLAW 
components are represented by the dotted (2008 data) and dashed 
(2010 data) lines. 



posure times of these observations were reduced with the 
duration of the eclipses contained in the data (500 s if a full 
eclipse was present, but less if only part of an eclipse was 
expected). Furthermore, Swift observations obtained within 
a 2-day time span were summed to improve the data statis- 
tics^ This seems justified, since the Chandra data do not 
reveal any spectral changes on such time scales (see Sec- 
tion [13}. 



2.4 Spectral models 

We fitted the spectra l data in the .5-10 keV energy range 
using XSPEC (v. 12.0; lArnaudlll996l ). This software package 
facilitates fitting a spectral model simultaneously to mul- 
tiple data files, which each have their own response and 
background files. As is common practise, we fit the XMM- 
Newton data with all spectral parameters tied between the 
different detectors (i.e., the model parameters are not al- 
lowed to vary independently between the PN and two MOS 
detectors). For all fits throughout this paper, we included 
the effect of neutral hydrogen absorption, Nu, along the 
line of sight using t he phabs model wit h the default XSPEC 
abundances (lAnders fc Grevessd 19891) and cross-sections 
|Balucinska-Church fc McCammorilll992r i. 

We first investigate the shape of the quiescent spectrum 
of EXO 0748-676 by considering the XMM-Newton observa- 
tion, which provides the highest statistics. A single absorbed 
powerlaw (powerlaw in xspec) provides an acceptable fit 
to the data (xt — 1.3 for 466 d.o.f.). However, the spectral 
index is unusually large for an X-ray binary (r = 4.7 ± 0.1) 
and suggests that the spectrum has a thermal shape. Using a 
simple absorbed blackbody model, bbodyrad, results in an 
adequate fit (xt = 1-2 for 466 d.o.f.), although the inferred 
emitting region has a much smaller radius than expected for 
a neutron star (~ 2 — 4 km for distances of 5 — 10 kpc). 



2 This is the case for obs IDs 31272037/38, 31272039/40 and 
31272043/44/45; see Table [T] 



Nevertheless, it is thought that radiative transfer effects 
in the neutron star atmosphere cause t he emergent spec- 
trum to deviate from a blackbody (e.g., IZavlin et al.lll996l ; 
iRutledge et al.. 1999). There are several neutron star atmo- 
sphere models av ailable within XSPE C , which yield equivalen t 
results (see e.g.. iHeinke et al.ll2006l ; IWebb fc Barretj|2007l 1. 
In the remainder of this work, we concentrate on fitting 
the data with a ne utron star atmosphere model nsatmos 
(|Heinke et al.ll2006l ). 

The NSATMOS model consists of five parameters, which 
are the neutron star mass and radius (Mns and -Rns), the 
effective temperature in the neutron star frame (i.e., non- 
redshifted; kT e fj), the source distance (D) and a normalisa- 
tion factor, which parametrizes the fraction of the surface 
that is radiating. We keep the latter fixed at 1 through- 
out this work, which corresponds to the entire neutron star 
surface emitting. The effective temperature as seen by an 
observer at infinity is given by kT£g = kT a g/(l + z), where 
1 + z = (1 — i? s /i?Ns)~ is the gravitational redshift factor, 
with _R S = 2GMns/c 2 being the Schwarzschild radius, G the 
gravitational constant and c the speed of light. 

The XMM-Newton data is well-fitted by an absorbed 
NSATMOS model (xt — 1-1 f° r 466 d.o.f.), although sig- 
nificant residuals above the model fit are present for en- 
ergies > 2 — 3 keV. We model this non-thermal emission 
by adding a powerlaw component, which significantly im- 
proves the fit (xt = 1.0 for 464 d.o.f.; an F-test suggests 
a ~ 1 x 10~ 14 probability of achieving this level of im- 
provement by chance). Chandra observations carried out in 
2008 mid-October, three weeks prior to this XMM-Newton 
observation, also indicated the presence of a non-thermal 
component in the qu iescent spectrum of EXO 0748-676 
IjDegenaar et al.ll2009l ). Whereas the Chandra data could 
not constrain the powerlaw index, the larger collective area 
of XMM-Newton provides better constraints for the fluxes 
under consideration. 

By using a combined NSATMOS and powerlaw model 
to fit the XMM-Newton data, we obtain a powerlaw index of 
r = 1.7 ± 0.5, i.e., in between the val ues of T = 1 and r = 2 
considered by iDegenaar et al.l (120091 ) . This fit furthermore 
yields JV H = (7±2) x 10 20 cm" 2 and Rns = 17.8±1 km, when 
fixing the neutron star mass to a canonical value of Mns = 
1.4 Mq and the distance to D = 7.4 kpc (the best est imate 
from type-I X-ray burst analysis; iGallowav et all20080 . The 
resulting powerlaw component contributes ~ 10 percent to 
the total unabsorbed 0.5-10 keV flux. This is lower than the 
~ 15 — 20 percent inferred fro m the Chandra observ ations 
performed in 2008 mid-October (JDegenaar et al.ll2009l ). The 
obtained hydrogen column density is consistent with values 
found for EXO 0748-676 d uring its outburst (N H - 7 x 
10 20 - 1.2 x 10 21 cm" 2 ; e.g.. ISidoli et~ai]|2005l ). 

The Chandra observations obtained in 2009 February 
and June are well-fitted by an absorbed NSATMOS model and 
do not require an additional powerlaw component. However, 
the 2010 April data shows evidence for such a hard tail, as 
significant residuals are present above the NSATMOS model 
fit for energies > 2 — 3 keV. If we include a powerlaw with 
photon index T = 1.7, as was found from fitting the XMM- 
Newton data (see above), this model component contributes 
~ 10, ~ 5 and ~ 15 percent to the total unabsorbed 0.5- 
10 keV flux for the data taken in 2009 February, June and 
2010 April, respectively. Fig. [3] compares the Chandra spec- 
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tral data obtained on 2008 October and 2010 April, showing 
that both spectral components decreased over the 18-month 
time span that separates the two observations. We found no 
spectral differences between the two separate exposures per- 
formed in 2009 February and therefore we tied all spectral 
parameters between these two spectra in the fits. 

The Swift data do not provide sufficient statistics to 
constrain the presence of a hard spectral component. We 
do include a powerlaw in the fits, but fix both the index 
and the normalisation of this component (see Section [3. 1|) . 
Since it is unclear how the powerlaw exactly evolves over 
time, we adjust the powerlaw normalisation for the Swift 
observations such that it always contributes 10 percent of 
the total unabsorbed 0.5-10 keV flux. After treating each 
Swift observation separately, we found that the thermal flux 
and neutron star temperature did not evolve significantly 
between consecutive observations. To improve the statistics, 
we therefore sum the Swift data into groups spanning ~ 1 — 4 
weeks of observations, resulting in exposure times of ~ 10 — 
20 ks (see Table [2]). The summed spectra were grouped to 
contain a minimum of 20 photons per bin. 



3 RESULTS 

3.1 Spectral fits 

As discussed in Section 12.41 the quiescent spectrum of 
EXO 0748-676 can be described by a combination of a neu- 
tron star atmosphere model and a non-thermal powerlaw 
tail. We fitted the Chandra and XMM-Newton data simul- 
taneously within XSPEC to a combined NSATMOS and power- 
law model subject to interstellar absorption, to explore the 
best-fit values for the neutron star mass and radius, source 
distance and hydrogen column density. We include the first 
set of Ch andra observations ob tained in 2008 October (dis- 
cussed in iDegenaar et al.l 120091 ) in the analysis. As before, 
we use the phabs model with the default XSPEC abundances 
and cross-sections to take into account the neutral hydro- 
gen absorption along the line of sight. The powerlaw index 
is fixed to V — 1.7 (the best fit- value obtained from XMM- 
Newton observations; see Section ^, 4[) . because there are not 
sufficient counts at higher energies in the Chandra spectra to 
allow this component to vary. The powerlaw normalisation 
is left as a free parameter. 

If the neutron star mass and radius are fixed to canon- 
ical values of Mns = 1.4 Mq and -Rns = 10 km, and in 
addition the source distance is fixed to D — 7.4 kpc, the 
hydrogen column density pegs at its lower limit (JVh = 0). 
When the distance is left to vary freely, the best-fit value 
is 4.6 ± 0.3 kpc, which is just outsi de the range obtaine d 
from X-ray burst analysis (5-8.3 kpc; lGallowav et alj |2008). 
Therefore, we choose to keep the distance fixed at 7.4 kpc, 
and instead allow the neutron star radius to vary. This way, 
we obtain best-fit values of Nu = (7 ± 1) x 10 20 cm -2 and 
R = 15.6 ± 0.8 km. If additionally the neutron star mass is 
left free to vary in the fit, this parameter is not strongly con- 
strained (Mns ~ 1-6 ± 0.6 Mo). In the final fits we choose 
to fix the neutron star mass to Mns = 1-4 M©, because 
otherwise the uncertainty in this quantity will dominate the 
errors of the other parameters. 

For the final spectral analysis, we fit all XMM-Newton, 
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Figure 4. Evolution of the bolometric flux (top) and effective 
temperature (bottom) of EXO 0748-676, deduced from Chan- 
dra /ACIS-S (black squares), Stui/t/XRT (grey triangles) and 
XMM-Newton/EPIC (black star) data. Multiple Swift obser- 
vations were summed to improve the data statistics (see Sec- 
tion!^. 



Chandra and Swift data with an absorbed NSATMOS plus 
powerlaw model, where Nh — 7 x 10 20 cm~ 2 , Mns = 
1.4 M , R NS = 15.6 km, D = 7.4 kpc and T = 1.7 are 
fixed, while the neutron star effective temperature is left as 
a free parameter. The powerlaw normalisation is left to vary 
freely for the Chandra and XMM-Newton observations, but 
fixed for the Swift data (so that this component contributes 
10 percent to the total unabsorbed 0.5-10 keV flux). We 
fit all data in the 0.5-10 keV energy range and deduce the 
absorbed and unabsorbed fluxes in this band. The thermal 
model fit is extrapolated to the energy range of 0.01-100 
keV to estimate the thermal bolometric flux. The results 
from fitting the X-ray spectra in this way are presented in 
Table [2] The effective temperatures and thermal bolometric 
fluxes derived from Chandra, Swift and XMM-Newton data 
are displayed in Fig. 2] Examination of Fig. U suggests that 
there is a small but discernible offset in the thermal flux 
and neutron star temperature as deduced from the different 
satellites. This is briefly discussed in Section f3. 31 



3.2 Lightcurve fits 

Fig. [4] clearly reveals a decaying trend in thermal flux and 
temperature. To investigate the decay shape, we fit the tem- 
perature curve with an exponential decay function of the 
form y(t) — a e ~^ t ~ t °'' T , where a is a normalisation con- 
stant, to is the start time of the cooling curve and r the 
e-folding time. Given the apparent offset between the differ- 
ent instruments (see Section [3. 3[) . we perform different fits 
to the Chandra and Swift data. We fix to to 2009 September 
5 (MJD 54714), which is in between the first non-detection 
by RX TE /PCA and the first Swift /XRT observation of the 
source (jDegenaar et al.l l2009V 



Table 2. Results from fitting the spectral data. 
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Satellite 



Date 



At 
(days) 



Pow. frac. 
(%) 



fcT cff 

(eV) 



pth 

bol 






(d.O.f.) 



Swift] 


2008-09-28 - 


- 2008-10-07 


1.9 


10 fix 


123.7 ±5.4 


1.31 ±0.22 


1.53 ±0.26 


10.0 ±1.7 


0.93 (8) 


Chandra] 


2008-10-12/13/15 


1.1 


20 ±3 


118.8 ±0.9 


1.23 ±0.02 


1.31 ±0.04 


8.6 ±0.3 


1.03 (175) 


Swift] 


2008-10-29 - 


- 2008-11-02 


2.2 


10 fix 


118.3 ±2.6 


1.10 ±0.09 


1.28 ±0.11 


8.4 ±0.7 


0.67 (14) 


XMM 


2008-11-06 




0.2 


7±2 


120.7 ±0.4 


1.14 ±0.01 


1.39 ±0.02 


9.1 ±0.1 


1.08 (467) 


Swift] 


2008-11-28 - 


- 2008-12-20 


11.0 


10 fix 


118.7 ±2.6 


1.11 ±0.10 


1.30 ±0.12 


8.5 ±0.8 


1.09 (14) 


Swift] 


2009-01-10 - 


- 2009-01-30 


9.8 


10 fix 


116.2 ±2.2 


0.99 ±0.07 


1.19 ±0.09 


7.8 ±0.6 


1.09 (18) 


Swift 


2009-02-13 - 


- 2009-02-23 


5.1 


10 fix 


117.2 ±2.4 


1.02 ±0.08 


1.23 ±0.10 


8.1 ±0.7 


0.50 (15) 


Chandra 


2009-02-23/25 


0.9 


12 ±4 


113.5 ±1.3 


0.91 ±0.03 


1.09 ±0.05 


7.1 ±0.3 


0.90 (139) 


Swift 


2009-03-01 - 


- 2009-03-16 


7.2 


10 fix 


115.6 ±2.3 


0.97 ±0.08 


1.17 ±0.09 


7.7 ±0.6 


0.79 (16) 


Swift 


2009-04-09 - 


- 2009-04-23 


7.1 


10 fix 


112.2 ±2.8 


0.86 ±0.09 


1.03 ±0.11 


6.8 ±0.7 


1.16 (10) 


Swift 


2009-05-07 - 


- 2009-06-05 


14.7 


10 fix 


114.2 ±2.3 


0.92 ±0.07 


1.11 ±0.09 


7.3 ±0.6 


0.88 (16) 


Chandra 


2009-06-10 




0.2 


4±3 


111.0 ±0.7 


0.75 ±0.01 


0.99 ± 0.03 


6.5 ±0.2 


1.19 (93) 


Swift 


2009-06-11 - 


- 2009-07-03 


11.5 


10 fix 


111.9 ±2.2 


0.75 ±0.07 


1.03 ±0.08 


6.7 ±0.5 


0.99 (17) 


Swift 


2009-07-18 - 


- 2009-07-31 


6.9 


10 fix 


110.5 ±2.1 


0.79 ± 0.06 


0.98 ± 0.07 


6.4 ±0.5 


0.42 (18) 


Swift 


2009-08-15 - 


- 2009-09-09 


12.7 


10 fix 


110.0 ±1.8 


0.78 ±0.05 


0.96 ± 0.06 


6.3 ±0.4 


1.30 (26) 


Swift 


2009-10-01 - 


- 2009-11-05 


17.4 


10 fix 


108.0 ±2.6 


0.69 ± 0.07 


0.88 ± 0.09 


5.8 ±0.6 


1.49 (11) 


Swift 


2009-12-21 - 


- 2010-10-01 


10.2 


10 fix 


109.4 ±2.0 


0.74 ± 0.06 


0.94 ±0.07 


6.1 ±0.5 


1.51 (19) 


Swift 


2010-02-12 - 


- 2010-03-13 


15.0 


10 fix 


109.4 ±2.0 


0.76 ± 0.06 


0.94 ±0.07 


6.1 ±0.5 


1.25 (20) 


Chandra 


2010-04-20 




0.2 


15 ±4 


108.6 ±1.1 


0.77 ±0.02 


0.91 ± 0.04 


6.0 ±0.2 


0.79 (91) 



Note. — The observations marked by a dagger were already discussed in iDegenaar et al.l 1 1200911 . but re- fitted in this work. These 
results were obtained by using a combined absorbed NSATMOS and POWERLAW model, where JVh = 7 x 10 20 cm" 2 , Afr-js = 1-4 Mq, 
^?NS = 15.6 km, D = 7.4 kpc and T = 1.7 were kept fixed. The quoted errors represent 90 percent confidence levels. Fx represents the 
0.5-10 keV total model flux and i^oi gives the 0.01-100 keV NSATMOS flux; both are unabsorbed and in units of 10 -12 erg cm -2 s . 
Lbol gives the 0.01-100 keV luminosity of the NSATMOS model component in units of 10 33 erg s _1 and assuming a source distance 
of D=7.4 kpc. At represents the time interval of the observations in days and the fractional powerlaw contribution is given in a 
percentage of the total unabsorbed 0.5-10 keV flux. 



The simple exponential decay, represented by the dotted 
lines in Fig.0 yields an e-folding time of 6121.7±2004.0 days 
for the Chandra data, but does not provide a good fit 
(xl = 6.0 for 2 d.o.f.). For the Swift lightcurve we find 
r = 5328.1 ± 674.7 days (xl = 0.5 for 12 d.o.f.). If we in- 



clude a constant offset (i.e., y(t) = a e 



-(t-t )/T 



+ b; solid 



lines in Fig. [5]), we obtain a better fit for the Chandra data, 
yielding a normalisation of a = 13.4 ± 0.2 eV, an e-folding 
decay time of r = 191.6 ± 9.7 days and a constant offset of 
b = 107.9 ± 0.2 eV (xl = 0.02 for 1 d.o.f.). For the Swift 
data we find a = 17.2 ± 1.8 eV, r = 265.6 ± 100.0 days 
and b = 106.2 ± 2.5 eV (xl = 0.34 for 11 d.o.f.), which is 
consistent with the Chandra fit. 

Although an exponential decay provides an adequate 
description of the dat a of EXO 0748-676, as has been found 
for ot her sources (e.g. JCackett et al.ll2006l : iFridriksson et al.l 
[2010J), mathematically a neutron star crust is expected 
to cool via a (broken) p owerlaw IjEichler fc Chend 1 19891 ; 
iBrown fc Cummind 2009). If we fit a single powerlaw of 
the form y(t) — A(t — to) B to the Chandra data, we find 
an index of B — —0.03 ± 0.01 and a normalisation of 
A = 134.4±1.0 eV (xl = 0.13 for 2 d.o.f.). For the Swift ob- 
servations we find B = -0.05 ±0.01 and A = 144.7 ±3.8 eV 
(xl = 0.4 for 12 d.o.f.). These powerlaw fits are indicated 
by the dashed lines in Fig. [5] 

A broken powerlaw also yields an acceptable fit to the 
Swift data (xl — 0.3 for 10 d.o.f.). We find a normalisa- 
tion of A = 135.0 ± 17.8 eV, a break at 166.0 ± 99.2 days 
and decay indices of —0.03 ± 0.03 and —0.06 ± 0.02 before 
and after the break, respectively. This fit is indicated by 
the dashed-dotted curve in Fig. [5] There are not sufficient 



Chandra observations to fit a broken powerlaw decay. We 
note that the shape of the decay curve of EXO 0748-676 
is not strongly affected by our choice of spectral parame- 
ters (Nh, Mns, -Rns, a nd Y) or assumed dis t ance (see also 
previ ous studies by e.g.. lWiinands et al.ll2004l ; ICackett et al.l 
120081) . 



3.3 Instrument cross-calibration 

The quiescent lightcurve presented in Fig. U shows indica- 
tions that the thermal flux and temperature inferred from 
the Chandra observations lie below the trend of the Swift 
data points. This possible shift (~ 6 percent for the flux 
lightcurve) may be due to cross-calibration issues between 
the two satellites. A study of the Crab nebula indeed re- 
vealed an offset between Chandra and Swift, whereas such a 
discrepancy was no t found between Swift and XMM-Newton 
(jKirsch et al.ll2005l ). This might be reflected in our results 
as well, since the XMM-Newton data point appears to line 
up with the trend indicated by the Swift data. However, 
our Chandra and Swift data points may also be (partly) off- 
set due to the fact that we cannot constrain the powerlaw 
component in the Swift data, which we therefore fixed to 
contribute 10 percent of the total 0.5-10 keV unabsorbed 
flux (see Section [2HI . 



4 DISCUSSION 

We discuss Chandra, Swift and XMM-Newton observations 
obtained after the cessation of the very long (~24 year) 
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Figure 5. Evolution of the effective temperature of EXO 0748-676 fitted to different decay functions (see Section I3.2II , The left image 
displays Chandra data and exponential decay fits both with and without a constant offset (solid and dotted line, respectively), as well 
as a decaying powerlaw (dashed curve). The right image shows Swift observations, where the dashed line is again a powerlaw fit, while 
the solid and dotted curves are exponential decays. In addition, this plot includes a fit to a broken powerlaw, which is represented by 
the dashed-dotted line. 



active period of EXO 0748-676. Fitting the spectral data 
with a neutron star atmosphere model NSATMOS, did not re- 
veal clear indications of a changing thermal spectrum during 
the fi rst five months of the quiescent phase (JDeeenaar et al] 
2009). However, now that the quiescent monitoring has ex- 
tended to 19 months (1.6 years), we find a significant de- 
crease in neutron star effective temperature from kT£§ ~ 
124 to 109 eV. The thermal bolometric flux was observed to 
decay from F b * ~ 1.5 x 10~ 12 to 0.9 x 10~ 12 erg cm" 2 s" 1 . 

In addition to a soft, thermal component, the Chandra 
and XMM-Newton observations show evidence for a hard 
powerlaw tail with index V — 1.7. The fractional contribu- 
tion of the hard spectral component to the total unabsorbed 
0.5-10 keV flux initially decreased from ~ 20 percent in 
2008 October to ~ 4 percent in 2009 June. However, obser- 
vations carried out in 2010 April suggest that the powerlaw 
fraction increased again to ~ 15 percent. Similar behaviour 
has been observed for severa l other quiesc ent neutron star 
systems (| Jonker et al] 120041 : Ijonkerl Bool), although oth - 
ers show more irregular behaviour i Fridriksson et al.ll2010f ). 
In Cen X-4, the powerlaw tail in the quiescent spectrum 
shows variations that appear to be linked to changes in the 
thermal component, possibly caused by low-level accretion 
jCackett et al.ll20ld ). 

The gradual decrease in thermal flux and neutron star 
temperature observed for EXO 0748-676 can be interpreted 
as the neutron star crust cooling down in quiescence after 
it has been heated during its long accretion outburst. Fig. [6] 
compares our data of EXO 0748-676 with the crust cooling 
curves observed for the neutron star X-ray binaries KS 1731- 
260, MXB 1659-29 and XTE J1701-462. This plot shows 
that the amount of cooling following the end of the outburst 
is markedly smaller for EXO 0748-676 than for the other 
three sources. We have observed our target over the first 



19 months after the cessation of the outburst and during this 
time the thermal bolometric flux has decreased by a factor of 
~ 1.7. In a similar time span, the thermal bolometric fluxes 
of KS 1731-260, MXB 1659-29 and XTE J1701-462 had 
decreased by a fact o r of ~ 3.5, 6 and 2.5, r espectively (see 
ICackett et all 120061 ; IFridriksson et~ai1 l201fj| ). The effective 
neutron star temperature of EXO 0748-676 has decreased 
by about 10 percent, compared to ~ 30, 40 and 20 percent 
for KS 1731-260, MXB 1659-29 and XTE J1701-462. 

Although the observed fractional changes in neutron 
star temperature and thermal bolometric flux are smaller 
for EXO 0748-676 than for the other three sources, the 
decay itself may not be markedly different. The quiescent 
lightcurves of KS 1731-260, MXB 1659-29 and XTE J1701- 
462 can be fit with an exponential decay function level- 
ling off to a constant value, yielding e-folding times of 
305 ± 50 , ~ 465 ± 25 and ~ 120 ± 25 da ys, respectively 



jCackett et alj|2008i : IFridriksson et al]|2010l ). For the Chan- 
dra data of EXO 0748-676, we find an e-folding time of 
~ 192 ± 10 days (see Section [3~2)l . These decay times provide 
a measure of the thermal relaxation time of the neutron star 
crust, which depends on the composition and structure of 
the lattice, the distrib ution of heating sources and the thick- 



ness of the crust (e.g ., Lattimer et al] 1994; Rut ledge et al] 



2002; Shtcrnin et al 



20071 ; iBrown fc Cumminell2009l 



iRutledge etHI (|2002T i and IShternin et "all |2007h cal- 
culate theoretical cooling curves for KS 1731-260, assum- 
ing different physics for the crust and core. These authors 
present simulations for both an amorphous crust and an 
ordered crystalline lattice. For the latter, the spread of 
nuclide charge numbers (Z) in the crust matter is small, 
which is referred to as a low level of impurities and re- 
sults in a highly conductive crust. A large number of im- 
purities gives an amorphous structure, which affects the 
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thermal properties of the crust and results in a low con- 
ductivity. In addition, iRutledge et al.l (|2002T ) explore stan- 
dard (i.e., slow) and enhanced neutrino cooling mech- 
anisms, yielding different core temperatures. Comparing 
our results on EXO 0748-676 with the decay shapes re- 
sulting from those calculations suggests that the neutron 
star has a highly conductive crust, simil ar to what has 
been i nferre d for the other t hree s o urces jW ijnand s et al.l 
20021. |2004 ICackett et~ai] l200d: IShternin et all 120071 ; 



Brown fc Cumminell2009l ; iFridriksson et alj|2010h . The fact 



that the decay curve of EXO 0748-676 is rather shallow 
may be explained in terms of a relatively small temperature 
gradient and thus lower thermal flux across the core-crust 
boundary (cf. the model curves for a highly c onductive crust 
and different core temperatures presented bv lRutledge et al.l 
120021 ) . This can be due to a combination of a warm neutron 
star core and a relatively low mass-accretion rate during 
outburst. 

The exponential decay fit to the Chandra data of 
EXO 0748-676 indicates that the neutron star crust might 
already be close to restoring equilibrium with the core. The 
fit results in a quiescent base level of 107.9 ± 0.2 eV, while 
we found a temperature of 108.6 ± 1.1 eV for the obser- 
vation performed in 2010 April. Prior to its last outburst, 
EXO 0748-676 was observed in quiescence with the EIN- 
STEIN observatory, displaying a 0.5- 10 keV u nabsorbcd 
flux of 8.4±£? x 10" 13 erg cm -2 s _1 l|Garcia fc Callananl 
1999). Our Chandra observations of 2010 April detected 
EXO 0748-676 at a 0.5-10 keV unabsorbed flux of (7.7 ± 
0.2) x 10 -13 erg cm* 2 s" 1 (see Table[2|. Assuming that the 
EINSTEIN detection caught EXO 0748-676 at its quiescent 
base level, this supports the idea that the crust has nearly 
cooled down. This would imply that t he neutron star cor e 
in EXO 0748-676 is relatively hot (cf. iHeinke et alj|2009h . 
suggesting that either standard cooling mechanisms are op- 
erating and that the neutron star is not very massive, or 
that the time-averaged mass-accretion rate of the system is 
very high due to a short recurrence time (see below). 

The energy deposited during outburst i s given by 
L nuc ~ (M)Qnuc/m u (e.g.. iBrown et all [19981 ; IColpi et all 
l200lf ). Here, Q nU c ~ 2 MeV is the nucle a r energy deposited 
per a ccreted baryon IjGupta et al.l 120071 ; lHaensel fc Zdunikl 
120081 ), m u is the atomic mass unit and (M) is the time- 
averaged accretion rate of the system. The latter can be 
expressed as (M) = (Mob) X iob/ircc, where (M b) is the 
average accretion rate during outburst episodes, t b is the 
outburst duration and tree is the system's recurrence time. 
The factor t h/t rcc represents the duty cycle of the system. 
The neutron star core is expected to be in a steady state, 
in which the energy radiated during quiescence balances the 
heat deposited during outburst. We can thus obtain an es- 
timate of the duty cycle of EXO 0748-676 by equating the 
heating and cooling rates. 

A neutron star cools primarily via photon radiation 
from the surface and neutrino emissions from the stellar 
core. If the lightcurve of EXO 0748-676 has indeed (nearly) 
levelled off, the bolometric luminosity emitted as photons is 
thus L-y ~ 6x 10 33 (D/7.4 kpc) 2 erg s _1 (as measured during 
the Chandra observation of 2010 April). The rate of neutrino 
emissions depends on the temperature of the neutron star 
core, which can be estimated from the effective surface tem- 
perature once the crust has thermally relaxed. A quiescent 



base level of kT^ ~ 108 eV (as suggested by exponential 
decay fits to the Chandra data), implies an effective surface 
temperature in the neutron star frame of fcT ff ~ 140 eV 
(~ 1.6 x 10 6 K), for a canonical values of M N s = 1-4 M Q 
and _Rns = 10 km (i.e., 1 + z = 1.3). Using the relation be- 
tween the effective surfa ce temperature and the inter ior tem- 
perature calculated by IBrown fc Cummind l|2009l ). yields 
Tcore ~ 1.3 x 10 8 K. For such a core temperature, the min- 
imum energy escaping the neutron star as neutrino's (i.e., 
assuming standar d core cooling) is L„ ~ ^q34-35 er g g -i 
jPage et 1 alj|2006h . 

Equating the energy losses via photon radiation from 
the neutron star surface (L-y) and neutrino emissions from 
the stellar core (L„) with the energy gained via crustal re- 
actions during outburst (L nuc ), suggests that EXO 0748- 
676 must have a time-averaged mass-accretion rate of 
(M) > 8 x 10 15 g s" 1 . During outburst, EXO 0748- 
676 displayed an average b olometric luminos i ty of ~ 6 x 
10 36 (D/7.4 kpc) 2 ergs" 1 (|Sidoli et al.l 120051 ; iBoirin et all 
|2007| ). Assuming that the accretion luminosity is given by 
Lace = (GMj<)s/RNs){Mob), this translates into a mass- 
accretion rate during outburst of (M b) ~3x 10 16 g s* 1 for 
a canonical neutron star with M = 1.4 Mq and R — 10 km(j 

If the crust has indeed thermally relaxed, the above es- 
timates show that EXO 0748-676 must have a duty cycle 
of > 30 percent to explain the observed quiescent bolomet- 
ric luminosity of ~ 6 x 10 33 (D/7.4 kpc) 2 erg s _1 in terms 
of thermal emission from the cooling neutron star (i.e., op- 
posed to continued accretion). The outburst of EXO 0748- 
676 started between 1980 May and 1984 July and the system 
returned to quiescence in 2008 September, i.e., t b = 24 — 28 
years. If the observed outburst is typical for the long-term 
behaviour of this source, the expected recurrence time is 
thus < 100 years. In case the neutron star cools via more 
efficient core neutrino emission processes, the recurrence 
time required to explain the observed quiescent luminosity is 
shorter (i.e., the duty cycle is higher). Although the above 
calculation is only a crude approximation (e.g., there is a 
significant uncertainty in the relation between the surface- 
and interior temperature of the neutron star, depending on 
the atmosph eric composition and the d epth of the light el- 
ement layer; IBrown fc Cummina 12009). it illustrates that 
EXO 0748-676 must have a high duty cycle if the cooling 
cur ve has indeed re ached its quiescent base lev el. 

IBrown et all l| 19981 ). IRutledge etafl (|200Ch and 
IColpi et al.l J200ll ') have suggested that EXO 0748-676 
continues to accrete in quiescence, because the quiescent 
luminosity inferred from the 1980 EINSTEIN observation 
is higher than predicted by standard cooling models. 
However, these conclusions are based on an assumed duty 
cycle of ~ 1 percent, but we have no a priori knowledge 
about this. Although we cannot exclude that the system is 
indeed accreting in quiescence, the above estimates show 
that a duty cycle of > 30 percent can explain the observed 
quiescent level of EXO 0748-676 as being due to thermal 



3 We note that EXO 0748-676 is an eclipsing system and there- 
fore part of the central X-ray flux may be intercepted from our 
line of sight. However, the X-ray burst behaviour of the source is 
consistent with th e mass-accretion ra te inferred from the observed 
X-ray luminosity iBoirin et al.l i2007). 
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emission from the cooling neutron star. A duty cycle of > 30 
percent is high, al though not unprec e dented for neutron star 
trans ients (e.g., IChen et al.l 1 19971 : iDegenaar fc Wiinandsl 
l2009h 

Recently, iBrown fc Cummina (J200y) demonstrated 
that the cooling of a neutron star crust is expected to fol- 
low a broken powerlaw decay. A break is predicted to occur 
due to a transition in the crystal structure of the crust mat- 
ter, and the slope before the break reflects the heat flux from 
the outer crustal layers. Therefore, we also fitted the neutron 
star temperatures obtained for EXO 0748-676 to a power- 
law and found decay indices of — 0.03±0.01 and — 0.05±0.01 
for the Chandra and Swift data sets, respectively. The Swift 
observations indicate that a possible break in the quiescent 
lightcurve may have occurred ~ 67 — 265 days after the ces- 
sation of the outburst (see Section \3. 2 1) . By fitting a broken 
powerlaw function, we obtain a decay index of —0.03 ± 0.03 
before the break, which steepens to —0.06 ± 0.02 thereafter. 
However, since these slopes are consistent with being equal, 
further observations are required to confirm whether a break 
has indeed occurred. 

The decay parameters that w e find for EXO 07 4 8-676 
are comparable to that obtained bv lFridriksson et al.l (| 20100 
for XTE J1701-462. These authors found that the quies- 
cent lightcurve breaks ~ 20 — 150 days post-outburst and 
report decay indices of ~ — 0.03 and ~ —0.07 befor e and 
after the break, respectively. iFridriksson et al.l (|2010l ) note 
that possible cross-calibration effects between Chandra and 
XMM-Newton might introduce small shifts that also allow 
a single powerlaw decay with slope ~ —0.05. The cooling 
curves of KS 1731-260 and MXB 1659-29 appear to have 
steeper de cays with indices o f ~ —0.12 and ~ —0.33, re- 
spectively (jCackett et al.ll2008l ). Due to the scarcity of data 
points it is unclear whether a break occurr ed in the qui- 
escent lightcurves of the latter two sources IjCackett et al.l 
|200S| ; IBrown fc Cumminell2009l ) . 

The powerlaw fits show no indications that the quies- 
cent lightcurve of EXO 0748-676 is levelling off. Thus, it is 
also possible that the neutron star temperature continues to 
decay further and that the core is cooler than suggested by 
the exponential decay fits and the 1980 EINSTEIN detec- 
tion. The relatively slow decrease of EXO 0748-676 might 
then reflect that the crust has a high conductivity, albeit 
lower than that of the neutron stars in KS 1731-260 and 
MXB 1659-29. Further observations are thus required to 
determine whether the neutron star crust in EXO 0748-676 
has nearly cooled down and to be able to draw firm conclu- 
sions on the crust and core properties. 
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Figure 6. Th e effective temperatu res of KS 1731-260 (green di- 
amo nds; from I Cackett et al.ll2006t) , MXB 1659-29 (red bullets; 
from ICackett et al.ll2006l. l200Sft ."xTE J1701-462 (grey crosses; 
from IFridriksson et alj|2010h and EXO 0748-676 (black squares, 
triangles and star). Exponential decay fits to the data of KS 1731— 
260, MXB 1659-29 and XTE J1701-462 are shown to guide the 
eye (green dashed, red dashed-dotted and grey dotted line, re- 
spectively). The two data points of XTE J 170 1-462 that lie above 
the decay fit are likely due to a temporary incr ease in the accre- 
tion r ate causing reheating of the neutron star iFridriksson et al.l 
|2010J) . 
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